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Abstract. This brief review provides an introduction to key ideas in the theory of 
incompressible magnetohydrodynamic (MHD) turbulence. 



1. Introduction 

Magnetohydrodynamic (MHD) turbulence has been studied in a num- 
ber of different parameter regimes, with important contributions from 
many authors, including [1-41]. This review focuses on incompressible 
turbulence in which the fluctuating velocities are comparable to or 
less than the Alfven speed [defined in equation (7)]. It is assumed 
that there is a mean magnetic field Bo that is at least as strong as 
the turbulent magnetic-field fluctuations. Sections 2 and 3 review the 
MHD equations and MHD waves, and section 4 describes wave-wave 
interactions. Section 5 reviews results on weak MHD turbulence, and 
section 6 treats strong MHD turbulence. Section 7 describes the role of 
dynamic alignment in the decay of MHD turbulence. 



2. Equations of incompressible MHD 

The governing equations of incompressible MHD are 

, (^ + v • w) _ -v L + §£) + -Lb ■ vb + (i) 



dB 

~dt 



V x (v x B) + ??V 2 B, (2) 

V • B = 0, (3) 

V • v = 0, (4) 



and 



p = constant, (5) 

where B is the magnetic field, v is the velocity, p is the thermal pressure, 
p is the density, and v is the viscosity. In the limit rj — ► 0, equation (2) 
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implies that field lines are like threads that are frozen to the plasma 
and advected at the plasma velocity v. 
Upon defining 

B = B z + 5B, (6) 



where Bqz is the uniform mean field, 

B 



v A z + <5b, (7) 



where va = Bo/y/4irp is the Alfven speed, the Elsasser variables 

a ± = v ± <5b, (8) 

and the total pressure 

n-£ + t (9) 

P 2 

and upon neglecting the dissipation terms (y = rj = 0), equations (1) 
and (2) can be rewritten as 



<9a + <9a^ 
- v A - 



dt dz 

and 



= -Vn-a _ Va + , (10) 



<9a da „ , N 

li - + «A ir = -Vn-a + -Va- 11 

dt dz 

In the absence of dissipation, equations (1) through (5) possess three 
quadratic invariants: the energy (1/2) / (v 2 + b 2 ) d?x, the cross helicity 
(1/2) / v • 5hd 3 x = (1/8) / (a + • a + — a - • a.~)d 3 x, and the magnetic 
helicity / A-B d 3 x, where A is the magnetic vector potential. The mag- 
netic helicity measures linkages of magnetic field lines and is associated 
with an inverse cascade [3]. The cross helicity describes the excess of 
one fluctuation type (a + or a - ) over the other, and is associated with 
the phenomenon of dynamic alignment (section 7). 



3. Waves and wave packets in incompressible MHD 

If v = (5b (i.e., a - = 0), the solution to equation (10) for V x a + is a 
wave that travels at speed va in the — z direction. Similarly, if v = — <5b 
(i.e., a + = 0), the solution to equation (11) for V x a~ is a wave that 
travels at speed va in the +z direction. These wave solutions are valid 
for arbitrarily large wave amplitude (when either a + = or a - = 0). 

There are two types of waves in incompressible MHD, the Alfven 
wave and the slow wave, with polarizations illustrated in figure 1. Both 
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ALFVEN-WAVE POLARIZATION 




Figure 1. For an Alfven wave, 5v and 5B are _L to both Bo and the wave vector k. 
For a slow wave, <5v and <5B are _L to k, but in the plane of k and B () . If k is nearly 
1 to Bo, then 5v and 8h are nearly along Bo for a slow wave. 



wave packet profile ! 



along a single field line: 1 1 1 1 I 1 1 1 



^B 



perturbed field line: 



if v = - 5b, then wave ^-r^ *~ B 



packet moves to the right: phase velocity = v A 

a - wave packet 
Figure 2. The profile of a ID Alfven wave packet. 

waves propagate along (or anti-parallel to) Bo at speed v\. A ID a - 
Alfven- wave-packet (with v = — <5b) is illustrated in figure 2. A cubical 
portion of a 3D wave packet is illustrated in figure 3. 

In general, a wave packet is associated with a displacement Ar of a 
field line, as in figures 2 and 3. For a small-amplitude wave packet that 
perturbs the field lines only slightly, 

Ar * r B±( ^ = ii<!^), (12) 

-£>o J-oo £>0 

where is the component of <5B in the xy-plane, and B^ is the 
Fourier transform of Bj_ in z. Thus, the field-line displacement is given 
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Figure 3. An example of a cubical portion of a 3D wave packet. 



by the k z = component of a wave packet. [16] If an Alfven wave 
packet's length along the field is An, its width across the field is Aj_, 
and its amplitude is B\ ± , then equation (12) implies that 

Aii-Bx, An cia i 
Ar~^i~^i, (13) 
Bo v A 

where a\ is the wave-packet amplitude in Elsasser variables. 



4. Wave-wave interactions in incompressible MHD 

For the nonlinear a 1 * 1 • Va^ terms in equations (10) and (11) to be 
nonzero, both a + and a fluctuations must be present at the same point 
in space. Wave-wave interactions can thus be thought of as collisions 
between oppositely directed wave packets [2]. When a + and a - wave 
packets collide, to lowest order in wave amplitude each wave packet 
follows the field lines associated with the other wave packet [11, 29]. 
For example, before a collision the value of a~ remains constant at 
a point that moves along Bo at speed v\. During the collision, the 
value of a - remains approximately constant at a point that moves at 
speed ~ along the field lines formed from the sum of Bo and the field 
fluctuations of the a + wave packet, with a pressure-induced correction 
to ensure V • a~ = [29]. This is illustrated in figure 4 for cubical 
portions of two colliding wave packets. In both cubes, <5B is upwards in 
the near-face of the cube, and downward in the rear-face of the cube, 
as in figure 3. 

The different displacements of different parts of a wave packet during 
a collision shear the wave packet in the plane perpendicular to Bo- For 
colliding wave packets of length An along Bq and width Aj_ across Bq, 
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BEFORE COLLISION: 



a - wave packet 



phase 
velocity 




field lines 

\ 



B, 



a+ wave packet 



AFTER COLLISION: wave packets have passed through each 

other and have been sheared 




B, 



phase 
velocity 



a - wave packet 



phase 
velocity 



Figure 4- During a wave-packet collision, the a (a + ) packet follows the field lines 
of the a + (a - ) packet. 



the change in the value of a 1 * 1 induced by the collision is roughly the 
duration of a collision, \\\/v\, times the magnitude of the nonlinear 



a^-Va 1 * 1 term. If la" 1 



~ a 



a\ , and if Aj_ < Am, then the fractional 



change x i n a induced by the collision is given by [13] 



X 



or, from equation (13), 



X 



Ar 



(14) 



(15) 



When % <C 1, the right-hand side of, for example, the a - wave 
packet in figure 4 is displaced in almost the same way as the left- 
hand side, since both sides encounter essentially the same a + wave 
packet, since the a + packet is only slightly altered during the collision. 
Changes to the profile of a wave packet along the magnetic field are 
thus weaker than changes in the profile of a wave packet in the plane 
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_L to Bo [8, 17, 18]. As a result, the cascade of energy to small Ay 
is much less efficient than the cascade of energy to small Aj_ in weak 
incompressible MHD turbulence, as discussed in the following section. 



5. Weak incompressible MHD turbulence 

In weak incompressible MHD turbulence, x C 1 and nonlinear wave- 
wave interactions are a small perturbation to a wave's linear oscillatory 
behavior. In this limit, the wave kinetic equation for the power spectra 
can be derived using weak turbulence theory [25]. When only Alfven 
waves are excited, the wave kinetic equation is of the form [25] 



dE ± (k) 



= J d 3 fcid 3 fe 2 [JV k)klik2 £^(k)-fl £)kl)ka £^(ki)]£; : F(k 2 ) 



dt 

x%k - <-<Wk-k!-k 2 ), (16) 
where k is the Fourier-space wave number, 

E ± (k)6(k + ki) = <a±(k) • a ± (k 1 )) ) (17) 

a 1 * 1 is the Fourier transform of a 1 * 1 , (...) indicates an ensemble average, 
the kernels A^ k kl k2 and i 3 k,k 1 ,k 2 are independent of E^ and given 
by [25], and 

w k = Tk z VA (18) 

is the frequency of small-amplitude a^-waves, from equations (10) 
and (11). Equation (16) contains only the leading-order terms in the 
weak-turbulence expansion, called the three-wave-interaction terms. 
The delta functions in the integrand give the three-wave resonance 
conditions, 

k = k!+k 2 , (19) 

and 

The integrand in equation (16) contains the product E^E^ and no 
terms of the form E ± E ± , consistent with the fact that only oppositely 
directed wave packets interact. Upon dividing equation (20) by v a, one 
finds that 

-k z = -k lz + k 2z (21) 
Combining equation (21) with the z-component of equation (19) yields 

k 2z = 0, and k z = k lz . (22) 
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If turbulence is not excited initially at some value of k z , three- wave 
interactions will not excite turbulence at that value of k z . Indeed, three- 
wave interactions do not transfer energy to larger values of k z , only to 

larger values of k± = \Jk% + ky [8, 25]. Also, the rate of energy transfer 

for is in essence proportional to the k z = part of . These 
conclusions remain valid when slow waves are also excited. [25] 

What physical meaning can be attributed to the k z = part of the 
power spectrum? It is simpler to answer this question by working with 
a ID spectrum (M, defined below) obtained by integrating E + over k x 
and k y . Assuming homogeneous turbulence, let 

(B + (x,y,0)-B + (x,y,z)) = g(z), (23) 

where B + is the fluctuating magnetic field associated with all a + wave 
packets (B + = a. + y/Wp), and let 

/oo 
e- lk * z g(z)dz. (24) 
-oo 

Equations (17) and (24) imply that M{k z ) = [p/32vr 4 ] / dk x dk y E + (k). 
Assuming that all wave packets have coherence length Ay along the 
magnetic field, coherence length Aj_ across the field, and rms fluctuating 
field B\ ± , one finds that 

M(0) = (/l B+ ( x >y>°)- B+ ( x >y> z ) dz )~ B L x \\- ( 25 ) 

If the wave packets have the Alfven-wave polarization, (5B _L Bo, then 
from equation (13), 

M(0) ~ B x± B Ar, (26) 

where Ar ~ B\ ± \\\/Bq is the field- line displacement associated with 
a single wave packet. Thus, for Alfven waves, M(0) [and similarly 
E^{k z = 0)] is a measure of the field-line displacements caused by wave 
packets. The dependence of the cascade rate on E Zf (k z = 0) indicates 
that energy cascade arises from the turbulent wandering of field lines, 
as described in section 4. [16-18, 30, 41] As illustrated by equation (26), 
the k z = part of the spectrum should not be equated with wave pack- 
ets that have infinite coherence length along the magnetic field [32, 41]; 
structures with finite Am lead to nonzero Ar, M(0), and E{k z = 0). 

The power spectrum of weak incompressible MHD turbulence can 
be estimated as follows [17, 18, 25, 28]. The fractional change \ m a 
wave packet of width A^ and length Am during one collision is <C 1. 
Since the changes induced by successive collisions add randomly, x 2 
collisions are required for an order-unity fractional change. The time 
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required for \ 2 collisions is \ 2 \/ v A, or ~ v A X\/ (a\^\\\). This is 
roughly the time T casca( je 

required for eddies of width Aj_ to pass their 
energy to eddies of width ~ Aj_/2. The cascade power e is roughly the 
energy per unit mass in wave packets of width Aj_, i.e. a\ /2, divided 

by ''"cascade • 

4 An 

e~^V- (27) 

Assuming that local interactions in wave-number space dominate the 
cascade, e is independent of Aj_ when Aj_ <C I and A^ S> d, where I is 
the outer scale (stirring scale) and d is the dissipation scale. Since Am 
is effectively constant, 

ax ± oc X 1 / 2 (28) 
[17, 18, 25]. If the one-dimensional power spectrum E\£,{k^_) is defined 



i j 



so that the energy of modes with k± in the interval (k±, 2ki), i.e. 
is given by J^ 1 E 1D (k ± )dk ± ~ k 1 E 1D (k 1 ), then 

E 1D (k ± ) oc kj 2 (29) 

[17, 18, 25]. 

Inserting equation (28) into equation (14) yields [18] 

X oc Xl 1/2 . (30) 

Thus, if the dissipation scale d is sufficiently small, then at sufficiently 
small Aj_ the value of x increases to 1, and the turbulence becomes 
strong. 



6. Strong incompressible MHD turbulence 

In strong MHD turbulence, for which % > 1, fluctuations are not waves, 
since they are significantly distorted by nonlinear interactions during 
a single wave period. Nor can propagation of fluctuations along the 
magnetic field be ignored, since, for example, for x ~ 1 propagation is in 
some sense equally important as wave- wave interactions. The combina- 
tion of propagation and nonlinear interaction determines the length of 
eddies in strong incompressible MHD turbulence as follows. [13, 29, 30] 
Consider a volume of width A^ across the magnetic field and length A 
along the field, with Ay ^ Aj_. The a~ fluctuations within this volume 
will be loosely described as a "wave packet" with the label Wl. It 
is useful to consider the "large-scale" magnetic field that is obtained 
by filtering out all magnetic fluctuations on length scales comparable 
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to or smaller than Wl. If one assumes that the dominant nonlinear 
interactions are between structures of similar size (local in fc-space), 
one can neglect the shearing of Wl by the large-scale magnetic field. 
Within a box that is large compared to Wl but small compared to the 
outer scale I, one can thus approximate the large-scale magnetic field as 
a uniform mean field denoted Bi oca i, as in figure 5. Prom equation (11), 
Wl propagates along Bi oca i at speed B\ oca x/ ^JA-Kp, while at the same 
time undergoing nonlinear interactions. Let W2 be a "wave packet" 
of a + fluctuations, also of width Aj_ and length Am, propagating an- 
tiparallel to Bi oca i while undergoing nonlinear interactions. Let a\ ± be 
the amplitude of a~ within Wl and also of a + within W2. Let Ay jCrit 
be the value of Ay such that % = 1 for Wl and W2: Ay iCrit = \±VA/a\ ± . 




Figure 5. In strong MHD turbulence, eddies are elongated along the local magnetic 
field. 

Let us suppose that Wl and W2 collide, and that Ay 3> Ay )Crit (i.e. 
X > 1), so that initially the a and a + fluctuations are coherent over a 
distance S> Ay crit along Bi oca i. After Wl and W2 have propagated into 
each other a distance Ay crit , as depicted in the bottom part of figure 6, 
the a - (a + ) fluctuations at the leading edge of Wl (W2) above point 
Q (P) have been altered by a factor of order unity by nonlinear interac- 
tions. As time proceeds, the a - fluctuations above points P and Q will 
evolve in different ways, because they are interacting with significantly 
different a + fluctuations: the a - fluctuations above point Q interact 
with the pristine, undistorted version of W2; the a - fluctuations above 
point P interact with a + fluctuations in W2 that have already been 
changed by the collision. As a result, the collision introduces structure 
with parallel coherence length Ay jCrit into the two wave packets. 
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According to these heuristic arguments, if % ls initially 3> 1 (but 
finite) , the combination of nonlinear interaction and propagation along 
field lines reduces x until x ~ 1. On the other hand, weak turbulence 
theory shows shown that if \ is initially <C 1, nonlinear interactions 
act to increase x until x ~ 1, as discussed in section 5. Thus, 1 is in 
some sense the "stable equilibrium value" of x m incompressible MHD 
turbulence. [13] 

When x ~ 1, the turbulence is strong. A theory for the spectrum and 
anisotropy of strong incompressible MHD turbulence was developed 
by [13]. The time 

^cascade f° r eddies of width A^ and length Aii to pass 
their energy to smaller eddies is ~ Am/ua ~ X±/a\ ± . Assuming local 
interactions in fc-space, the cascade power e is given by a^ ± /T casca d e ~ 
a| ± /Aj_. Assuming that e is independent of Aj_, 

oa ± oc A^ 3 , (31) 

which corresponds to 

E 1D (k ± ) oc kj_ 5/3 . (32) 
The condition % ~ 1 combined with equation (31) yields [13] 

Ay oc X 2 / 3 . (33) 
If the turbulence is stirred at scale I, and if a/ ~ v\, then [13] 

Ay = ZV3 A 2/3_ (34) 

Equations (31) through (34) apply to both Alfven-mode eddies and 
slow-mode eddies [13]. However, Alfven-mode fluctuations dominate 
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the shearing of both Alfven-mode eddies and slow-mode eddies [13]. 
The reason is that the nonlinear a 1 * 1 • Va T terms are much larger for 
Alfven modes than for slow modes, since fluctuations vary most rapidly 
across the magnetic field, and a 1 * 1 is perpendicular to the local magnetic 
field for Alfven modes, but nearly along the magnetic field for slow 
modes when Aj_ <C An. 

Equations (31) through (34) are consistent with the numerical re- 
sults of [27, 36]. The numerical simulations of [26] agree with equa- 
tion (32), but do not address local anisotropy. On the other hand, the 
numerical simulations of [29] are consistent with equation (34) , but find 

E±D(k±) oc kj 3 ^ 2 . Thus, although equations (31) through (34) are sup- 
ported by physical arguments and some direct numerical simulations, 
areas of uncertainty remain. 

7. Decay and dynamic alignment 

If only one type of fluctuation, either a + or a , is present, nonlinear 
interactions vanish. Such a state is called "maximally aligned," since 
either v = b or v = — b. Decaying "turbulence" that is maximally 
aligned decays only on a viscous or resistive time scale, and can thus 
be long-lived. 

If decaying turbulence initially contains both a + and a~ fluctua- 
tions, but an excess of one over the other, then it decays to a maximally 
aligned state, in which only the initially predominating fluctuation type 
remains. [4, 6, 9, 25, 29, 36, 41] The reason for this is the following. 
Whether the turbulence is weak or strong, isotropic or anisotropic, the 
energy-cascade time of a 1 * 1 fluctuations at a perpendicular scale Aj_ is 
proportional to some power of the energy in a^ fluctuations at scale A j_ . 
If turbulence is excited with an excess of a + waves, the decay time for 
the a - waves is short compared to the decay time of the a + waves. The 
a waves thus decay faster than the a + fluctuations, with a discrepancy 
in decay rates that increases in time. It has been suggested by [37] 
that this mechanism also operates in compressible MHD turbulence, in 
which case it provides a potential mechanism for long-lived turbulence 
in astrophysical environments. 
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